Recombining WMAP: 
constraints on ionizing and resonance radiation at recombination 
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We place new constraints on sources of ionizing and resonance radiation at the epoch of the 
recombination process using the recent CMB temperature and polarization spectra coming from 
WMAP. We find that non-standard recombination scenarios are still consistent with the current 
data. In light of this we study the impact that such models can have on the determination of 
several cosmological parameters. In particular, the constraints on curvature and baryon density 
appear to be weakly affected by a modified recombination scheme. However, it may affect the 
current WMAP constraints on inflationary parameters like the spectral index, n s , and its running. 
Physically motivated models, like those based on primordial black hole or super heavy dark matter 
decay, are able to provide a good fit to the current data. Future observations in both temperature 
and polarization will be needed to more stringently test these models. 



I. INTRODUCTION 

The recent measurements of the Cosmic Microwave 
Background (CMB) flux provided by the Wilkinson 
Microwave Anisotropy Probe (WMAP) mission pj 
have truly marked the beginning of the era of precision 
cosmology. In particular, the position and amplitude 
of the detected oscillations in the angular power spec- 
trum of the CMB are in spectacular agreement with 
the expectations of the standard model of structure 
formation, based on primordial adiabatic and nearly 
scale invariant perturbations 2] . Assuming this model 
of structure formation, an indirect but accurate mea- 
surement of several cosmological parameters has been 
reported in agreement with those previously indi- 
cated (see e.g. and |f|). 

However, beside a full confirmation of the previous 
scenario (with a sensible reduction of the error bars!) 
the WMAP data set is also hinting towards a modi- 
fication of the standard picture in several aspects 0- 
In particular, the low CMB quadrupole JaJ3] possible 
scale dependence of the s pec tral index y| Q , high op- 
tical depth (see e.g. |ToL llll IT^l possible deviations 
from flatness 01 an d dark energy (|l4|) have already 
produced a wide interest. Further, the relatively high 
X 2 per degree of freedom of the fiducial model suggests 
systematics and/or new physics could be considered. 

It is therefore timely, with the increased precision of 
the CMB dataset to further test the standard scenario 
and to investigate possible deviations. 

Here we investigate possible deviations in the mech- 
anism to which CMB anisotropics are higly depen- 
dent: the process of recombination. 

The recombination process can be modified in sev- 
eral ways. For example, one could use a model inde- 
pendent, phenomenological approach such as in 
where models are specified by the position and width 
of the recombination surface in redshift space. Here 
we will instead focus on theoretically motivate mecha- 
nisms based on extra sources of ionizing and resonance 



radiation at recombination (see e.g. 1.6] )■ While 
the method we adopt will be general enough to cover 
most of the models of this kind, we remind the reader 
that there are several other ways in which to modify 
recombination, like, for instance, by having a time- 
varying fine-structure constant ([13). In our analy- 
sis, we will not cover the constraints on these models 
from WMAP as they have been recently investigated 
in(0)._ 

Following the seminal papers hj, \m detailing 
the recombination process, further refinements to the 
standard scenario were developed [2l| allowing pre- 
dictions at the accuracy level found in data from 
the WMAP satellite and the future Planck satellite 
p2l I23I ] . With this level of accuracy it becomes con- 
ceivable that deviations from standard recombination 
maybe be detectable [lll2ll25|. 

The letter proceeds as follows: in section [H] we re- 
view models which can produce deviations from the 
standard recombination scenario. In IIIII we describe 
how these deviations might affect the CMB temper- 
ature and polarization power spectra and conduct a 
likelihood analysis using the recent CMB data from 
WMAP. In particular, we will study the impact that 
a modified recombination scheme can have on sev- 
eral cosmological parameters. In II VI wc draw together 
the implications of the analysis, placing constraints on 
current theories of recombination. 



II. BEYOND STANDARD RECOMBINATION 

In the standard recombination model [h! l2pj | the 
net recombination rate is given by 



dx e 
~~dt 



std 



c 



a c nx 2 e — b c (l — x e ) exp 



AB 



(1) 



where x e is the ionization fraction, a c and b c are the 
effective recombination and photo-ionization rates for 
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principle quantum numbers > 2, and AB is the differ- 
ence in binding energy between the I s * and 2 nd energy 
levels. In addition to the single Ly-a transition with 
wavelength X a , there is also 2 photon decay from the 
meta-stable 2s level, with decay rate Ai S 2 S . The con- 
tribution of this process is reflected in the parameter 
C, 



c = 



1 + KA ls2s n ls 
1 + K(A ls2s + b c )m 



K = 



Al 



8nH(z) 



(2) 



The standard hydrogen recombination scenario can 
be simply extended in two ways with the addition of 
Ly-a and ionizing photons [l6|. 

The extra contributions can be related to the 
baryon number density through the efficiency func- 
tions e a and Si, respectively, 



dn a . , . . dn^ 

— = e a (z)H(z)n, — 



e t {z)H(z)n (3) 



where H(z) and n are the Hubble expansion parame- 
ter and mean baryon density (neutral H and protons) 
respectively. We leave photons produced by two body 
interactions, for example of annihilation of WIMPS, 
with drii/dt — uwiMPnwiMP f° r future investiga- 
tion. 

The addition of extra Ly-a and ionization photons 
in equation |J3J| adjusts the recombination rate from 
the standard model in equation Q 



dx e 
~~dt 



dx e 
~dt 



\std 



CeiH - (1 - C)e a H. (4) 



There are a range of physical mechanisms which 
could generate additional photons, we review them 
briefly here, see and references therein for a fuller 
discussion and derivations. 

In general e a and £j are functions of redshift, de- 
pendent upon the particle decay model that is being 
employed, and the size of the decay lifetime, t x in 
comparison to the lifetime of the universe t^. 

Primordial black hole (PBH) decay 26] at the mo- 
ment of recombination can be described by a model 
of the form 



jC 2 exp -C" 



C 



1 + z 

1 + Z de , 



(5) 



where j can be a or i and Zd ec is the redshift of pho- 
ton decoupling. The authors of [25| use c a « 0.3 and 
Ci pa 0.13 however these values are dependent on as- 
sumptions about the PBH mass and number density 
distributions, which in turn rely on the nature of the 
inflationary spectrum [23] and the possibility of accre- 
tion [28|. 

Electromagnetic cascades from general particle de- 
cay can been described by the parameterization [25j ] 
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FIG. 1: Evolution of characterizing the number density 
of extra ionizing photons added during the recombination 
epoch (bottom panel) and the evolution of the associated 
ionization fraction x e (top panel) for the WMAP best fit, 
fiducial case (full line, Ei = 0) and for sample models, de- 
scribed in the text: PBH decay with coefficients as in pfll ] 
(dot-long dash), topological defect particle release (dot- 
short dash), short lived SHDM decay with B = 10 5 (short 
dash), and = 10 6 (long dash) corresponding to z(t x ) ~ 4 
and 5 respectively. The ionization fraction for the topo- 
logical defect model lies almost on top of the fiducial case. 
The parameter measuring additional resonance photons, 
e a , is of the same order ase^ in each of these models. 



with r = 1/2 for particle release from topological de- 
fects, r = — 1 for decay of super heavy dark matter 
(SHDM), and where the normalization is the EGRET 
energy density [2ij, although one could imagine us- 
ing a code such as DARKSUSY [3(j to get a more 
precise normalization value. The proportionality con- 
stant Olz^)] reflects an alteration required for de- 
cay time dependence in the normalization of SHDM 
scenarios. For short lived SHDM with t x < to, 



e[ z {t a 



and equals 1 for all 



- exp \il + z[t x ]f' 2 

other cases. Short lived SHDM can therefore have 
substantially higher values for e ita . 

Note that in both particle and BH decay the models 
predict ~ e a . We demonstrate the size and evolu- 
tion of Ci for scenarios described above in Fig. ^ 

The combined effect on reionization from mecha- 
nisms such as those described above would obviously 
require a sum over contributions from all sources. Our 
aim is to find an upper limit on the overall contribu- 
tion, independent of source, and as such we adopt a 
simple parameterization using constant, effective val- 
ues for Ei and e a . 

In Figs. [21 and |3| we plot the temperature (TT) 
and temperature-polarization cross correlation (TE) 
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FIG. 2: CMB TT spectra for various values of e a (top) and 
£i (bottom) showing the increasing suppression and shift of 
the peaks as one increases the number of extra resonance 
and ionizing photons, respectively, in comparison to the 
WMAP best fit fiducial model with e l = e a = (full line). 
Spectra are normalized to Cs7- 



FIG. 3: CMB TE spectra for various values of e a (top) 
and Si (bottom) against a fiducial case with e a — £i = 
(full line). In each case a reionization redshift of z ri — 7 
is used for comparison. 



power spectra for several values of e a and e^. Qual- 
itatively, increasing e a broadens the epoch of recom- 
bination and lowers the redshift when the optical 
depth drops through unity. Delaying recombination 
increases the angular diameter distance at last scat- 
tering and therefore shifts the first peak to lower I. 
The delay also increases the optical depth, suppress- 
ing the height of the peaks in comparison to the low I 
plateau and decreasing the ratio of the second to first 
peaks. 

Increasing £j introduces a plateau in ionization frac- 
tion, preventing a trend towards full recombination. 
Boosting the number density of ionizing photons, in- 
creases the optical depth significantly more than a 
similar increase in Ly-a photons, resulting in a pro- 
nounced suppression on the first peak height. Subse- 
quently we would expect much tighter constraints on 
£i from the TT and TE spectra. 

It is interesting to note that introducing a large 
extra ionizing component can give a large cross- 
correlation on large scales, similar to that generated 
by early reionization. However it is unable to account 
for both the TT and TE observations simultaneously 
because of the peak suppression in the TT spectrum. 



III. LIKELIHOOD ANALYSIS 

Our analysis method is based on the computation 
of a likelihood distribution over a grid of precomputcd 
theoretical models. 



We restrict our analysis to a flat, adiabatic, A-CDM 
model template computed with a modified version 
of CMBFAST 1^3, sampling the parameters as fol- 
lows: Slcdmh 2 = u>cdm = 0.05, ...0.25, in steps of 0.01; 
n b h 2 = io b = 0.009, ...,0.030, in steps of 0.001 and 
h = 0.55, 0.85, in steps of 0.05. The value of the 
cosmological constant A is determined by the flatness 
condition. Our choice of the above parameters is mo- 
tivated by Big Bang Nucleosynthesis bounds on w& 
(both from D j^l and 4 He + 7 Li 33]), from super- 
nov ae |34i | and galaxy clustering observations (see e.g. 
[351 136| ). From the grid above we only consider mod- 
els with age of the universe to > 11 Gyrs. We vary 
the spectral index of the primordial density perturba- 
tions within the range n s = 0.8, ..,1.2, we allow for 
a possible (instantaneous) reionization of the inter- 
galactic medium by varying the reionization redshift 
5 < z r i < 25 and we allow a free re-scaling of the 
fluctuation amplitude by a pre-factor of the order of 
C 87 , in units of C^ ORM = 1.9fiK 2 Finally, we let 
e a and £j vary as follows: 10~ 4 < e a < 10 2 , and 
< 10~ 4 < E; L < 10 2 in logarithmic steps. 

The theoretical models are compared with the 
recent temperature and temperature-polarization 
WMAP data using the publicly available likelihood 
code pLOl l37|| . 

In Fig. 2] we plot the likelihood contours in the 
e a — Ei plane showing the 2a and 3<r contours. As we 
can see there is no strong correlation between the two 
parameters. 

Marginalizing over all the remaining nuisance pa- 
rameters we obtain the following constraints: e a < 
lO 5 and ei < 1Q- 1 - 2 at 95% c.L. 
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FIG. 4: Likelihood contour plot in the e a — £j plane show- 
ing the 2a and 3<r contours. 
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TABLE I: Difference in \ 2 from the best fit standard re- 
combination model for a range of non-standard recombi- 
nation scenarios discussed in section [H] using WMAP TT 
and TE data. For the SHDM models we consider a range 
of decay lifetimes, t x , measured in terms of the lifetime of 
the universe at a redshift, z(t x ). 



As we can see, despite the high precision of the 
WMAP data, substantial modifications to the recom- 
bination process might be present. In Table I we 
demonstrate this for models discussed in section [HJ 
Ionizing photons from PBH decay are still in good 
agreement with the data, to within la. Similarly, as- 
suming the normalization prescription given in [25j . 
long lived SHDM particles are not constrained well 
by the data, however there is a sharp cut off in the 
likelihood for particles that would decay on timescales 
shorter than t[z = 6). 

Since a modified recombination scheme can still be 
in agreement with the data, it is interesting to study 
the impact these modifications might have on the con- 
straints of several parameters determined in the stan- 
dard analysis. In Fig. we demonstrate the ben- 
efits of having both the TT and TE spectra avail- 
able in breaking potential degeneracies with spatial 
curvature. With pre- WMAP uncertainties in the TT 
spectrum peak heights relative to the plateau, a 'de- 
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FIG. 5: Figures demonstrating how a potential degeneracy 
between flat and spatially curved scenarios with non-zero 
e a is broken with WMAP TT and TE spectra. A WMAP 
best fit fiducial model with z r i = 17 and Qk = (full line) 
is shown against a model with Qk = 0.6 and e a = 10 3 
and z r i = 17 (dashed line) that prior to the WMAP data 
release had some degeneracy with the standard scenario. 
In particular note how the addition of the TE spectrum 
strongly breaks the degeneracy, since the open model has 
a strongly suppressed reionization peak. 



generacy' between curvature and e a existed A 
non-zero value of e a can shift the TT spectrum peak 
positions in an open model so that they have a simi- 
lar a fiducial A CDM case. However the accuracy of 
WMAP data over the first and second peaks, now al- 
lows allows such an open model to be distinguished by 
the peak suppression produced by increasing e. More- 
over, the reionization peak is strongly suppressed in 
the open scenario so the two scenarios are further dif- 
ferentiated by the inclusion of TE data. 

In Fig. El we probe degeneracies that are less easily 
distinguished; projecting the constraints on e a and £$ 
in function of the baryon density, the spectral index 
and its running dn s /dlnk at fco = O.ObhMpc^ 1 . As we 
can see, the constraints on uji, are weakly affected by 
our simple modifications to the recombination process. 
The value of uj b = 0.023 ± 0.001 is unaffected by the 
inclusion of e a or This can be explained by the 
high precision measurements of Sachs- Wolfe plateau 
and the first 2 acoustic peaks made by the WMAP 
satellite, which breaks the degeneracy between wt, and 
Ei found in Seager et al. Increasing e a or £j, however, 
can have an important impact on the determination of 
n s . Namely, a modified recombination will have effects 
similar to those of an early reionization process, low- 
ering the small scale anisotropics and allowing greater 
values ofn s to be in agreement with the data. 
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FIG. 6: Correlations between e a and e% with several cos- 
mological parameters. The redshift of reionization is fixed 
at z r i — 12 and the Hubble parameter is h — 0.68. One 
can see that while the baryon density is robust to alter- 
ations in e a and e%, both the scalar spectral index and its 
running are sensitive. 



Finally, it is interesting to study the correlation 
with a possible running of the spectral index. In 
Peiris et al. (fH), a ~ 2a evidence for a negative 
running of the spectral index has been found. This is 
not expected in most common inflationary scenarios 
and several mechanisms have been proposed to ex- 
plain the effect (see e.g. 9j). However, as we can 
see from the plots at the bottom of Fig|SJ a delayed 
recombination will in general not solve this problem, 
but it will enhance it towards a more negative run- 
ning (in the plots we assume n s = 0.93). A possible 
solution might come from a speeding up of recombina- 
tion rather than delay, through having concentrations 
of baryons in high density, low mass clouds, for ex- 
ample |24|. The accelerated recombination, because 
of the resultant premature reduction in the ionization 
fraction, can be parameterized by an effective e a < 
with £j > still. As we see in FigjJJ negative values of 
e a are more consistent with a zero running, although 
arc likely not sufficient to fully explain the observed 
phenomenon. 




Log(-eJ 

FIG. 7: Correlation between a negative e a and a running 
of the spectral index dn/dlnk at k = 0.05hMpc . The 
redshift of reionization is fixed at z ri — 12, the Hubble 
parameter is h — 0.68, n s — 0.93 and e% — 



IV. CONCLUSIONS 

We have probed the upper bounds that can be 
placed on the contribution of extra Ly-a and ioniz- 
ing photon-producing sources through their effect on 
recombination and subsequently on the CMB TT and 
TE spectra. We have found that, adopting a simple 
parametrization using constant, effective values for e a 
and Si, the WMAP data contraints e Q < 10 and 
e t < HT 1 - 2 at the 95% level. 

We have found that the WMAP data is able to 
break many of the worrying degeneracies between e a , 
Si and the more standard cosmological parameters. In 
particular, the constraints on curvature and baryon 
density appear to be weakly affected by a modified 
recombination scheme. However, it may affect the 
current WMAP constraints on inflationary parame- 
ters like the spectral index n s and its running. 

Physically motivated models, like those based on 
primordial black hole or super heavy dark matter de- 
cay, are able to provide a good fit to the current data. 
Future observations in both temperature and polar- 
ization, from the next WMAP release and the Planck 
satellite p5|. will be needed if we are to more strin- 
gently test these models. 
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